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ABSTRACT 

We used the Wide-field Infrared Camera on the Canada-France-Hawaii telescope to observe four 
transits of the super-Earth planet GJ 1214b in the near-infrared. For each transit we observed in two 
bands nearly-simultaneously by rapidly switching the WIRCam filter wheel back and forth for the 
duration of the observations. By combining all our J-band (~1.25 fim) observations we find a transit 
depth, analogous to the planet-to-star radius ratio squared, in this band of (i?p,//i? Hl ) 2 =1.338±0.013% 
- a value consistent with the optical transit depth reported by Charbonneau and collaborators. How- 
ever, our best-fit combined Ks-band (~2.15 /jm) transit depth is deeper: (i?pif s /i?») 2 =1.438±0.019%. 
Formally our Ks-band transits are deeper than the J-band transits observed simultaneously by a fac- 
tor of (RpKs/ Rp j) 2 — I -072±0. 018 - a 4cr discrepancy. The most straightforward explanation for our 
deeper Ks-band transit depth is a spectral absorption feature from the limb of the atmosphere of the 
planet; for the spectral absorption feature to be this prominent the atmosphere of GJ 1214b must 
have a large scale height and a low mean molecular weight. That is, its atmosphere would have to be 
hydrogen/helium dominated and this planet would be better described as a mini-Neptune. However, 
recently published observations from 0.78 - 1.0 /im, by Bean and collaborators, show a lack of spectral 
features and transit depths consistent with those obtained by Charbonneau and collaborators. The 
most likely atmospheric composition for GJ 1214b that arises from combining all these observations 
is less clear; if the atmosphere of GJ 1214b is hydrogen/helium dominated then it must have either a 
haze layer that is obscuring transit depth differences at shorter wavelengths, or significantly different 
spectral features than current models predict. Our observations disfavour a water-world composition, 
but such a composition will remain a possibility for GJ 1214b, until observations reconfirm our deeper 
Ks-band transit depth or detect features at other wavelengths. 

Subject headings: planetary systems . stars: individual: GJ 1214 . techniques: photometric- transits 
- infrared: planetary systems 



1. INTRODUCTION 

Astronomers have been waiting for sometime for 
a planet remotely similar to our own Earth that 
could be readily investigated with current instruments. 
Such an object was recently announced with the sem- 
inal discovery of the sup er-Earth planet GJ 1214b 
(jCharbonneau et al. I I2009f) with the MEarth telescope 
netw ork (jNutzman fc Charbonneau 1 120081 : llrwin et al. I 
2009). Although not the first transiting super-Earth 
announced - CoRoT-7b arguably holds that honour 
(|Leger et al. I 120091 : iQueloz et all l2009t iPont et al. I 
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l20T(it) - GJ 1214b is in many ways more interesting as 
it offers the opportunity for advantageous follow-up to 
constrain its planetary characteristics. With a mass of 
6.55 M ffi and a radius of 2.68 R m , GJ 121 4b's density 
(p~1.87 g cm -3 ; iCharbonneau et all 120091 ) is less than 
that of the terrestrial planets of our solar system and 
therefore GJ 1214b may have a significant gaseous at- 
mosphere. Also, as it transits a low mass star, its equi- 
librium temperature (T e9 ~500K assuming a low Bond 
albedo) is much more hospitable than CoRoT-7b, and it 
has a much more favourable planet-to-star radius ratio; 
as a result, if there are significant spectral features in its 
atmosphere then they should be detectable with current 
instruments. 

As there are no super-Earth analogues in our solar 
system, it is a pressing question whether the burgeon- 
ing c lass of planets with minimum masses b elow 10 
(e.g. lUdrv et al. I [20071 IMavor et al. Il2009albl ) are pre- 
dominantly scaled-down Neptunes, with large helium- 
hydrogen atmospheric envelopes, or scaled-up terrestrial 
planets with atmospheres predominately composed of 
heavier molecules. Fortunately, GJ 1214b is an ideal can- 
didate to answer such questions: iMiller- Ricci fc Fortnev I 
(2010) showed that due to GJ 1214b's advantageous 
scale-height and planet-to-star radius ratio it should have 
readily observable water and methane spectral features 
across the infrared spectrum if its atmosphere is com- 
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posed primarily of hydrogen and helium. Conversely, 
if its atmosphere is composed predominantly of heav- 
ier molecules, then the resulting smaller scale height will 
mute the spectral features and current instruments will 
return transit depths consistent with the depths mea- 
sured in the optical bv ICharbonneau et al. I ()2009h . 

Furthermore, recent theoretical work has shown that 
measurements to constrain the composition of the 
gaseous atmosphere of GJ 1214b will also constrain the 
plan et's bulk composition. | Rqgers &: Seagerl ([2010ft and 
later iNettelmann et al. I (12010ft showed thatthe observed 
mass and radius ( Charbo nneau et al. I 12009ft can be 
equally well fit by either a significant rocky core/mantle 
and a hydrogen-rich atmosphere, or a "water- world" with 
a small water-rich core and a significant steam atmo- 
sphere. Searches for spectral features will differentiate 
between these two scenarios; detections of prominent 
spectral features will argue for the former scenario of 
a rocky core to go along with the hydrogen/helium-rich 
gaseous pla netary envelope. 

Recently, iBean et al. I (|2010f ) have performed just such 
transit spectroscopy o bservations u s ing th e FORS2 in- 
strument on the VLT; IBean et al~l (12010ft obtained 11 
spectrophotometric light curves from 0.78 - 1.0 [im that 
show consistent t ransit depths with one anoth er. By 
comparing to the IMiller-Ricci h Fortnevl (|2010l ) atmo- 
spheric models they were able to show that the lack of 
observed spectral features suggest that GJ 1214b must 
either have a high mean molecular weight and is likely 
a water- world, or its atmosphere is hydrogen/helium 
dominated with hazes or clouds high in the atmosphere 
that obscure the expected spectral features shortward 
of — 1 micron. An additional possibility could certainly 
be that GJ 1214b's atmosphere is more complicated 
than expected, and its atmosphere could still be hydro- 
gen/heli um dominated with diff e rent spectral features 
than the IMiller-Ricci fc Fortnevl (|2010ft models suggest. 

We have also performed broadband transmission spec- 
troscopy observations searching for GJ 1214b's spectral 
features from —1-2.5 microns using the Wide-field In- 
frared Camera (WIRCam) on the Canada-France-Hawaii 
telescope (CFHT). We have already successfully demon- 
strated the precision of WIRCam on CFHT in the near- 
infrared through our detections of the secondary eclipses 
and ther mal emission for TrE S-2b and TrES-3b in the 
Ks-band (|Croll et al. I l2010alTbh. and for WASP- 12b in 
the J, H & Ks- bands (jCroll et all 120 lift . Here we re- 
port observations of several transits of the super-Earth 
GJ 1214b in three bands with WIRCam on CFHT; 
for each transit we observed near-simultaneously in two 
bands to allow for accurate comparisons of the transit 
depths between these two bands. We observe an in- 
creased transit depth in the Ks-band as compared to the 
J-band depth, likely indicative of absorption near —2.15 
fim. The only way to achieve an absorption feature this 
prominent is if GJ 1214b has a large scale height, a low 
mean molecular weight, and thus its atmosphere is hy- 
drogen/helium dominated. We discuss below the likely 
possibilities for the atmospheric m ake-up of GJ 1214 b 
that result from a combinat ion of the lBean et al. I (|2010f ) . 
ICharbonneau et al~l (|2009ft and our own data. Our re- 
sults disfavour a water-world composition, but such a 
composition is possible if our Ks-band point is simply 
an outlier; such a composition will remain a distinct 



possibility until further observations either confirm our 
increased Ks-band depth or detect spectral features at 
other wavelengths. The observations to date are most 
qualitatively consistent with a hydrogen/helium domi- 
nated atmosphere that is either hazy or one with more 
complicated spectral features than our current models 
suggest, such as an atmosphere where non-equilibrium 
chemistry plays a significant role. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the transit of GJ 1214b on four occasions 
with WIRCam on CFHT. For each transit we observed 
GJ 1214 in two bands nearly simultaneously by rapidly 
switching the WIRCam filter wheel back and forth be- 
tween the two filters. We observed GJ 1214 alternat- 
ing between the J (—1.25 /irn) and Ks (—2.15 \mi) fil- 
ters on the evenings of 2010 June 27, 2010 August 15 
and 2010 September 22. On the evening of 2010 Au- 
gust 7 we observed GJ 1214 alternating between the J 
and CH 4 On filters (-1.69 ptm); these 2010 August 7 ob- 
servations were of a partial transit only, as the airmass 
of GJ 1214 increased during transit and during egress 
reached an airmass of —2.3, close to the hard limit of the 
telescope, preventing further observations. For the 2010 
June 27 transit the airmass was low for the duration of 
the observations, while for the 2010 August 15 and 2010 
September 22 transits the airmass increased throughout 
the observations and exceeded an airmass of two by the 
end of the out-of-transit baseline. For these latter two ob- 
servations we noticed that the data quality significantly 
degraded as the airmass rose above 2.0; for the analysis 
that follows we therefore exclude the data in the out-of- 
transit baseline with an airmass greater than two for the 
2010 August 15 and 2010 September 22 o bservations. 

We observed GJ 1214 in Staring Mode (jDevost et al. I 
l20Toh where we use the full WIRCam array with its 
21'x21' field-of-view and do not dither for the duration of 
our observations. The only exception to our normal star- 
ing mode practice was the aforementioned filter change. 
The exposure times and defocus amounts for our various 
observations were: 3.5 s and 2.0 mm for Ks-band, 4 s 
and 1.2 mm for the CH40n-filter, and 3.5 s and 1.8 mm 
for J-band, respectively. The filter change during our 
observations induced an additional overhead of 20 sec- 
onds to switch the filter wheel, as well as an additional 
60 seconds to take two acquisition images to fine-tune 
the target position and reinitiate guiding. We observed 
in data-cubes of 12 exposures to increase the observing 
efficiency. We took 12 exposures (one guide-cube) in a 
single filter, before performing the filter change and ob- 
serving an additional 12 exposures in the other filter. 
The resulting duty cycles were: 22% for our J and Ks- 
band observations, and 19% for our CH40n and J-band 
observations. 

The data were reduced and aperture photometry was 
performed on our target star and numerous reference 
stars. We used apertures with radii of 17, 18 and 15 pix- 
els for our Ks, J and CH4O11 photometry, respectively; 
the associated inner and outer radii for the sky annuli 
were 21 and 29 pixels for our Ks-band and CH40n pho- 
tometry, and 22 and 30 pixels for our J-band photometry. 
We preprocess our data using the I'iwi pipeline, designed 
specifically for WIRCam imaging. We performed differ- 
ential photometry on the target star with between 13-25 
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Figure 1. The root-mean-square of our out-of-eclipse photometry 
for our various data-sets. The dashed line in each panel displays the 
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(solid line) after the subtraction of their respective background trends 
one over the square-root of the bin-size expectation for Gaussian noise. 



reference stars for our photometry in our various bands; 
for further details o n the method we refer the reader to 
ICroll et al. I (pOrO a.b). As can be seen in Figure [TJ the 
out-of-transit photometry after subtraction of the back- 
ground trend bins down near the Gaussian noise ex- 
pectation with increasing bin size in all cases. Our data 
is therefore not seriously affected by time-correlated red- 
noise. We set our errors for our transits equal to the 
root mean square of the data outside of transit after the 
subtraction of a linear trend; for our J and CH^On pho- 
tometry on 2010 August 7 we scale up the errors from 
the RMS by a factor of | as this data is of a partial 
transit only, and the egress of transit occurs at very high 
airmass, which we found to be correlated with degraded 
precision with our other data-sets. The resulting light 
curves for the various observations are displayed in the 
top panel of Figure [2J 

3. RESULTS 

We fit each of our data-sets with a IMandel fc Agol I 
(2002) transit model, with the depth of transit, 
(i? p /i?») 2 , and the best- fit mid-transit time as free pa- 
rameter^. For several of our GJ 1214 datasets we 
noticed obvious trends with time after our differential 

8 We quote the barycentric Julian date in the ter restrial time 
standard using the routines of (Eastman ct al. 2010). 



photometry was performed, similar to the trends no- 
ticed in several of ou r exis ting WIRCam datasets (e.g. 
ICroll et al. I l2010allE I2MTI ). We cannot rule out that 
these trends are intrinsic to the target star and could 
be due to, for instance, long-term stellar variability as a 
result of rotational modulation. However, the frequency 
with which we observe such trends suggests that most of 
these trends are likely systematic in nature. We there- 
fore refer to these trends as background trends, and we 
fit our datasets with linear or quadratic backgrounds of 
the form: 

B f = 1 + b x + b 2 dt + b 3 dt 2 (1) 

where dt is the time interval from the beginning of the ob- 
servations and bi, bi and 63 are fit parameters. To deter- 
mine whether a quadratic (63) term is justified to account 
for the background trend , we calculate the Bayesian in- 
formation criterion fBIC: iLiddle Il200"7l ). and ensure that 
the BIC is lower with the inclusion of the quadratic term. 
That is, the reduction in x 2 must be sufficient to justify 
the extra degree of freedom. Only the August 15 Ks- 
band data warranted a quadratic term (63). 

We employ Markov Chain Mont e Car lo (MCMC) 
fitting as described in ICroll I (|2006[ ) and ICroll et al. 1 
(2010a) using chains with 5xl0 6 steps. There are four 
free parameters for each data-set: (i? p /i?*) 2 , the best-fit 
mid-transit time, b\, and 62- We also fit our data with 
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Figure 2. CFHT/WIRCam photometry bracketing the transit of GJ 1214 obtained nearly simultaneously on the nights denoted in the 
figure in the J-band (left panels) and the Ks-band and the CH40n filter (right panels). The top set of panels (A) show the data and the 
best-fit transit model displayed with the appropriately coloured line. The second set of panels (B) show the residuals of the data and the 
models from a transit model with a depth given by the best-fit value of the individual J-bands transits at left. Thus the (B) left panel 
simply displays the residuals from the best-fit transit model from that night's data. The (B) right panel displays the residuals of the Ks 
and CFLiOn filter data and model from a transit with a depth equal to the best-fit J-band transit depth observed on that night nearly 
simultaneously and using the appropriate Ks or ClUOn limb-darkening coefficients. The third set of panels (C) are the same as the (B) 
panels, except with the data binned every twelve points. The bottom panels (D) display all the data obtained simultaneously in the J and 
Ks-bands from the (C) panels binned every twelve points [that is the top three transit curves in the ABC panels]; the orange dot-dash 
line displays the difference from the best-fit J-band transit model using the combined J and Ks-band best-fit depths (Q. The errors in 
this panel are calculated from the standard deviation of the points within each bin. The best-fit Ks-band transits display increased transit 
depths, while the best-fit CH40n transit displays a smaller transit depth than the J-band transits observed simultaneously. In the top 
three panels the different data-sets are offset vertically for clarity. Note each set of panels has a different vertical scale. 
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the "residual-permutation" method (W inn et al. 1 120091) 
where the residuals to the best-fit data are shifted and re- 
fit, thus preserving the time correlation of any red-noise 
in the data, so as to investigate the impact of time corre- 
lated systematics. We find similar results to our MCMC 
analysis, although in general the errors are slightly larger. 
The increased size of the errors otherwise is likely due 
to low number statistics as in this method for N data- 
points one can only generate 27V- 1 light curve permuta- 
tions (~200-350 iterations for each of our data-sets). For 
these reasons we quote our MCMC results henceforth. 

We employ a quadratic limb-darkening law and ob- 
tain o ur limb-darkening parameters from iClaret et al. I 
(1998) for the J, Ks and CH^On filters - we adopt their 
H-band limb-darke ning paramete r s for o ur CH^On filter 
observations. The IClaret et al. I (|1998[ ) limb-darkening 
parameters are calculated t hrough fits to the PHOEN IX 
stellar atmosphere models (jHauschildt et ainil997al |bl). 
We employ the parameters ci and C4 from non-linear 
limb darkening laws, as quoted in Table [TJ For the in- 
put values that we use to generate the limb-darkening 
parameters we use approximations to the measured stel- 
lar effective temperature {T e ff — 3026 K), and log of 
the stellar surface gravity (log g = 4.991 [CGS units], 
iCharbonneau et al. 112009 ), of T eff = 3000 K and log g 
= 5.0. We adopt the period and ephemeris of GJ 1214b 
given m IBean et al. I (f2010h . We use the inclination 
and semi-major axis to s tellar radius raticH determined 
from an analysis from the IBean et al. I (|2010f ) white-light 
curves (J. Bean electronic co mmunication). All other pa - 
rameters were fixed at the ICharbonneau et al. I ()2009f ) 
values. Our best-fit MCMC transit and background fits 
for each of our individual transit datasets are listed in the 
"Ind." rows of Tableland displayed in the appropriately 
coloured lines of Figure 

As discussed below (Q we note a small, but signifi- 
cant difference in the transit depths we measure in the 
J and Ks-bands. To confirm that this difference is sig- 
nificant, and not due to any uncertainty in the limb- 
darkening coefficient^ we also fit our three J & Ks-band 
transits simultaneously while fitting the two quadratic 
limb-darkening coefficients for each band (c2j, C4J, C2Ks, 
and CiKs)- We place a priori constraints on the limb- 
darkening coefficients; these priors on the limb-darkening 
coefficients are Gaussian with a s tandard devia t ion o f 
0.05 from the values derived from IClaret et al. I (1998) 
(and listed in the "Ind." rows of TableUJfor the J and Ks- 
bands). For each of our three transits we fit the Ks and 
J-band data-sets with their own background terms (e.g. 
bu, b2j, biKs, and b2Ks for the 2010 June 27 transit). 
We fit for the three mid-transit times of our three transits 
(to allow for possible transiting timing variations), but 
ensure that this value is held in common between the J 
and Ks-band transits observed simultaneously. We fit for 
the J-band transit depths, but fit the Ks-band transits 
with a depth that is a consistent factor [{Rpks/ Rpj) 2 } 
greater, or less than, the J-band transit observed nearly 
simultaneously for all three transits. By utilizing this 

9 An inclination of i=88.94° and a semi-major axis to stellar 
radius ratio of a/i? s =14.97. 

10 We also us e d qua dratic limb darkening coefficients from 
ILester fc Ncilson (20(53) an d a model without limb-darkening. In 
both cases our Ks-band transit depths were significantly deeper 
than the J-band depths. 



fitting method, in addition to the methods discussed be- 
low in £14.31 we can directly ascertain how much deeper 
the Ks-band transits are than the J-band transits. The 
advantage of fitting for the depth of each of the J-band 
transits individually, rather than fitting them with a con- 
sistent depth, is that this allows for small variations in 
the J-band depth that could arise from rotational modu- 
lation (as discussed in £|4.2p , while still directly fitting for 
the ratio of the Ks to J-band transit depths. We fit for 
24 parameters overall, and the best-fit results are listed 
in the rows marked "Joint" analysis in Table [TJ 

3.1. WIRCam non-linearity correction 

We also ensured that any difference in the transit depth 
from the J to Ks-bands was not due to an ineffective non- 
linearity correction. During the Fiwi preprocessing step, 
a non-linearity correction is applied to correct the count 
levels for pixels that approach saturation. Near satura- 
tion this non-linearity correction can be as large as 10%. 
At the maximum count levels recorded in a pixel of the 
aperture of our target star, GJ 1214, during our obser- 
vations, the detector is well below its saturation level, 
and the WIRCam detector is approximately 3-5% non- 
linear at these count levels. The vast majority of the 
pixels in our target star and reference stars apertures are 
illuminated to much lower levels and are expected to be 
non-linear at the 1-3% level. If this non-linearity correc- 
tion was applied ineffectively then this could cause a sys- 
tematic offset in our measured transit depths; although 
this discrepancy was expected to be much smaller than 
the difference in the transit depth from J to Ks that 
we measure, we nonetheless demonstrated this was the 
case by reprocessing and reanalyzing our 2010 Septem- 
ber 22 transit data in the J and Ks-bands without apply- 
ing the non-linearity correction. Any deviations in the 
pixel count values from the current non-linearity correc- 
tion, will be more than an order of magnitude smaller 
than the effect induced by not applying the non-linearity 
correction whatsoever. Not employing the non-linearity 
corrections, as expected, leads to shallower transit depths 
than when the non-linearity correction is applied. How- 
ever, the ratio of the transits depths from Ks to J are 
near identical whether the non-linearity correction is, or 
is not, applied. Overall, as this test should create a vari- 
ation much larger than one due to an ineffective non- 
linearity correction, it is safe to conclude that the greater 
Ks-band than J-band transit depth does not arise from 
the non-linearity correction. 

4. DISCUSSION 

4.1. GJ 121Jfh's transit depth in the near-infrared 

We display our best- fit transit depths in Figure [3] and 
Table [TJ The J-band transit depths are largely consis- 
tent with one another and are also consistent, or at most 
insignificantly shallower, t han t he depths re p orted by 
(ICharbonneau et al. 1 12009) and IBean et al. 1 (|2010f ) in 
the optical and very near-infrared. The Ks-band transits 
also display similar depths to one another. However, the 
Ks-band transits appear to be deeper than the J-band 
transits; this is a small effect, but is clearly visible in the 
bottom panel of Figure [5] where we present the residu- 
als of our observations from the best-fit J-band transit 
depths observed nearly simultaneously. 
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Table 1 

CFHT/WIRCam near-infrared transit depths of GJ 1214b 



Date 
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Fit 
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a Fit Type stands for the joint or individual (Ind.) analyses. 

The CH4O11 transit depth, on the other hand, appears 
to have a similar transit depth to the the J-band transit 
observed simultaneously on 2011 August 7. As this is a 
partial transit only, and as much of the transit and the 
egress of transit occurs at very high airmass, caution is 
warranted in any robust comparison of the J to CH4O11 
transit depth and to other wavelengths. 

By combining all the Ks-band and J-band transits, we 
find the weighted means and the associated errors on 
the transit depths are (i? pJ /i?*) 2 =1.338±0.013% for J- 
band, and (R pKs / i?*) 2 = 1.438±0.019% for Ks-band. As 
we only have one partial transit of GJ 1214 in the CELtOn 
filter, the depth in that band is simply the value from 
the 2010 Aug. 7 transit: (i^cH 4 On/-R*) 2 =l-290± u ; u g%. 
We determine the error on the weighted mean of our J- 
band and Ks-band points by determining the weighted 
error of all our observations in that particular band and 
then scaling that error upwards by a factor of To 
determine ( we calculate the x 2 of all our data in a sin- 
gle band compared to a model with a consistent transit 
depth equal to the weighted mean of the transit depths 
in that band; we then scale up the errors to ensure the 
reduced x 2 is equal to oncEH- The Ks-band data-points 
are consistent with one another, so only the J-band er- 
rors are scaled upwards. The resulting value is £j=1.02 
for the J-band photometry, so this suggests that both the 
Ks and J-band weighted errors are already appropriately 
sized, or close to it. Overall, this analysis suggests that 
our Ks-band and J-band transit depths are inconsistent 
with one another; the Ks-band transit depth is deeper 
than the J-band depth with 5c confidence. 



4.2. The effect of stellar spots on transit observations of 

GJ 1214b 

iCharbonneau et al. I (|2009| ) reported that GJ 1214 is 
an active star and displays longer term variability with 
a period of several weeks at the 2% level in t he MEarth 
bandpass. More recently, iBerta et al. I (|20100 presented 
and analyzed MEarth photometry of GJ 1214 from 2008 
to 2010 and observed long-term photometric variability 
at the 1% level with a period of approximately 50 days. 
This variability is presumably due to rotational modula- 
tion from spots rotating in and out of view. As the long- 
term p hotometric monitoring presented in IBerta et al. I 
( 2010) ends in 2010 July (in the midst of the observa- 
tions we present here) we assume the more conservative 
limit of 2% variability for our calculations henceforth. 

Transit observations obtained at different epochs may 
show small differences in the transit depth due to ro- 
tational modula tion arising from bo t h occulted and un- 
occulted spots (iCzesla et al. 1 120091 IBerta et all 120101: 
iCarter et al. II2011D . In the case of unocculted spots, if 
the 2% observed rotational modulation represents the 
full range from a spotted to unspotted photospher43 
then we may expect measurements of the transit depth 
of GJ 1214b will vary by as much as 0.03% of the stel- 
lar signal in the MEarth bandpass (assuming an unspot- 
ted transit depth of 1.35%) from observations taken at 
epochs spanning the maximum and the minimum of the 
observed rotational modulation. On the other hand, oc- 
culted spots will cause small brightening events during 
the transit that may lead one to underestimate the true 
transit depth. Thus, for transit depth measurements ob- 
tained at different epochs, such as our own, it is possi- 



11 lAnd rac et al. (201Cj) notes that there are several hidden as- 
sumptions one should be careful to address when applying reduced 
X 2 to one's data; we feel the method we apply here should be 
useful nonetheless as a first-order approximation to indicate the 
appropriate size of the weighted error bars 



There is no reason to expect we ever observe a hemisphere 
of the star free of spots altogether, and in deed the analysis of 
ICarter et al. I H2011I) and IBerta et al. I (120101 ') suggest we very well 
may not, which would lead to larger transit depth differences 
fr om epoch to epoch. Exact analytical expressions are available 
in ICarter et~aT~l |[_0T_ ). 
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Figure 3. WIRCam transit observations of the super-Earth GJ 1214b. The WIRCam response functions are displayed inverted at 
the top of the plot using the black dotted lines at an arbitrary scale; these are from left to right: J-band (~1.25 /im), the CH^On 
filter (~1.69 /ira), and Ks-band (~2.15 /ira). Our J and Ks-band observations are indicated by the red diamonds for our 2010 June 27 
observations, the blue diamonds for our 2010 August 15 observations, and the magenta diamonds for our 2010 September 22 observations. 
Our 2010 August 7 J-band and CH-jOn filter observations are indicated by the green diamonds. We offset our measured J and Ks 
observations slightly in wavelength for clarity. For the Ks and J-band we display th e weighted mean and error of t he observations in these 
bands with the black diamonds. The orange solid and cyan dashed curves are the [Miller- Ricci & Fortneyl 1120101 ) GJ 1214b atmospheric 
mode ls for solar metallicity h y droge n/helium dominated atmospheres (with and without met hane, respectively). Th e grey dotted curve 
is the [Miller- Ricci & Fortnex] 112010 ) models for an H^O/steam atmosphere. We integrate the Millcr-Ricci & Fortncy (2010) atmospheric 
models over the WIRCam response functions and display these values in the appropriately coloured solid triangles. 



ble that the variability caused by rotational modulation 
creates small differences in the measured depths for ob- 
servations taken at different epochs. 

However, for our data obtained nearly simultaneously 
in two different bands, the effect of spots should be re- 
duced for the following reasons. First of all, as the stellar 
rotation period of GJ 1214 appears to be mu ch longer 
(jCharbonneau et al. Il200l iBerta et al. 1120101 ) than the 
~1 hour duration of the transit, the spot pattern should 
be essentially static during a single transit. Secondly, 
even if the star is very spotted during our own observa- 
tions, the difference between the transit depths measured 
nearly simultaneously in our two bands will be minute. 
This small difference will arise due to the differing ra- 
tio of the Planck function of the spot and the star due 
to their different temperatures; however, this difference 
will be muted as we move into the near infrared. For 
instance, assuming that GJ 1214 has spots 500 K cooler 
- a value supposedly consistent with another M4.5 dwarf 
(jZboril et al. 1 l2003h - than GJ 1214's ~3000 K effec- 
tive temperature " ijCharbonneau et al. I [20091. the 2.0% 
variability due to rotational modulation in the MEarth 
bandpass (^780 nm), will translate into 1.5%, and 1.0% 
variability in J-band and Ks-band, respectively. Assum- 
ing the unspotted transit depth is 1.350% in these bands 
then the maximum transit depths from unocculted spots, 
which would results from measurements at the minimum 
flux of the observed rotational modulation, would be 
1.369% in J and 1.364% in Ks-band. The variation in 
the transit depth between the J and Ks-bands that we 
observe is both much larger than this predicted effect 
due to starspots, and also would serve to create a deeper 
transit in J-band, rather than Ks-band; we, of course, 
observe the opposite phenomenon. 

A deeper Ks-band than J-band transit could arise 
from spots along the transit chord that are oc- 
culted during the observations. Occultation of spots 
by a planet will create anomalous brightenings dur- 
ing transit, as has bee n observed for th e transit- 
ing planets HD 189733b (jPont et al. 1 120071 ). TrES-lb 



(|Rabus et a l. Il2009t [ Dittmann et aT~1 12010D. and more 
recent ly for GJ 121 4b (IBean et al . Il2010t IBerta et al. I 
120101: iCarter et al~l |2011j). However, due to the near- 
simultaneous nature of our photometry, it is unlikely 
that occulting spots could account for the variation in 
Ks to J-band transit depth that we observe, again due 
to the small difference in the Planck functions of the 
spot relative to the star between our two bands. For in- 
stance, occulting a spot 500 K cooler than the surround- 
ing photosphere that is 30% of the planetary size, will 
lead to transits that are 0.07% shallower for the dura- 
tion of the occultation than the presumed 1.35% transit 
depth in J-band. However the Ks-band transit depth will 
also be 0.05% shallower; the 0.02% relative difference be- 
tween the two bands expected from a spot occultation is 
much less than the observed transit depth difference we 
observe. Lastly, as spot occultations will lead to over- 
all shallower transit depths, one would require the true 
transit depth of GJ 1214b to be deeper than that ob- 
served in our own J-band o bservations and in the op- 
tical and very nea r-infrared (jCharbonneau et al. 1 120091 : 
IBean et al. II2010T ). For these reasons we find it unlikely 
that the transit depth difference we observe arises from 
occulted or unocculted spots. 

4.3. A larger transit depth in Ks-band than J-band 

Due to the aforementioned possible variations in the 
transit depths from epoch to epoch induced by spots, a 
more straightforward method to compare the depths of 
transits in our bands is to directly compare the depth of 
the transit in one band to the depth obtained simulta- 
neously in another - that is (Rpks/ Rpj) 2 for each one 
of our transits observed simultaneously in the Ks and J- 
bands. We ignore our data observed on 2010 August 7 
for this analysis, as GJ 1214 was observed in the J-band 
and the CH40n filter, rather than in J and Ks. 

We attempt to measure the fraction that the Ks- 
band transits are deeper than the J-band transits, 
{RpKs/ Rpj) 2 , by two methods. In the first method, 
we display the best-fit MCMC transit depth of our Ks- 
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Figure 4. The measured Ks-band transit depth divided by the 
measured J-band transit depth for the "Ind." analysis of our var- 
ious observations (see the legend). The horizontal dotted line de- 
notes the weighted mean of the Ks divided by J-band transit depth, 
{RpKs/Rpj) 2 • The solid horizontal solid line denotes the value if 
the Ks-band transits were the same depth as the J-band transits. 
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Figure 5. The measured Ks-band transit depth divided by the 
measured J-band transit depth from our "Joint" analysis. The ver- 
tical solid line denotes the best-fit value, while the dotted vertical 
lines indicate the 68% credible regions. 

band photometry divided by the J-band transit depth 
from our individual analysis ( "Ind." ) in Figure 2] The 
associated errors displayed in this figure are generated by 
propagating through the associated errors on the individ- 
ual best-fit MCMC transit depths as displayed in Table 
[TJ The weighted mean and error of these data indicates 
that the Ks-band transits are deeper than the J-band 
transits by a factor of: (R PKs / R PJ ) 2 =1.072±0M8. The 
associated errors do not need to be scaled up, as the 



reduced 



is near one for a comparison of our three 

2 



(RpKs/Rpj) 2 data-points as compared to the weighted 
mean of these observations; specifically x 2 =1.92, which 
is reasonable given the two degrees of freedorrP^I By an- 
alyzing the individual transits, our Ks-band photometry 
displays a deeper transit depth than our J-band photom- 
etry observed nearly simultaneously with a confidence in 
excess of 4a. 

In the second method we use the "Joint" MCMC anal- 
ysis described above where we simultaneously fit the 
three three data-sets that observe in Ks-band and J- 
band simultaneously. We fit each J-band transit with 

13 We again note there are several hidden assumptions one 
shou ld be aware of when applying reduced \ 2 as documented in 



an independent transit depths, but assume all the Ks- 
band transits were a consistent factor deeper (or shal- 
lower) than the J-band transits observed nearly simulta- 
neously, {RpKs/ Rpj) 2 ■ We also fit the limb-darkening 
parameters after applying a Gaussian a priori assump- 
tion as described in [J3] We display this ratio for our 
"Joint" analysis in Figure [SJ we measured this frac- 
tion as (Rpks/Rpj) 2 = 1.063t°;°2?- This distribution is 
not a perfect Gaussian, and therefore, according to this 
method our Ks-band transits are deeper than our J-band 
transits with greater than 3<r confidence. 

Both methods return reasonably similar results, and 
argue for a deeper Ks-band transit depth than J-band 
transit depth with a confidence in excess of 3<r. We quote 
the value of (Rpks/ Rpj) 2 from our "Ind." results hence- 
forth. 

4.4. WIRCam transit depths suggest a low mean 
molecular weight 

The transit depths from our CFHT/ WIRCam photom- 
etry suggests that GJ 1214b should have a large scale 
height, low mean molecular weight and thus a hydro- 
gen/helium dominated atmosphere. These conclusions 
arise from the fact that our Ks-band transit depths are 
deeper than our J-band transit depths by a factor of 
(i? PK;j /.Rp,/) 2 =1.072±0.018. This corresponds to a rel- 
ative change in the radius of GJ 1214b from the Ks-band 
to the J-band of 1.04%, or ~610 km compared to GJ 
1214b's radius of ^17070 km. Using its equilibrium tem- 
perature of 7^^560 K (the value obtained assuming zero 
bond albedo) this amounts to absorption increasing the 
radius of GJ 1214b by ~2 atmospheric scale heights for a 
hydrogen gas envelope (H2). As, the molecular weight of 
H2O is approximately nine times greater than hydrogen 
gas, a spectral feature this prominent assuming a wa- 
ter world composition would require an increase in the 
planetary radius of ~20 atmospheric scale heights. This 
suggests that an atmosphere dominated by light elements 
is much more probable than one composed of heavier el- 
ements. 

The change in transit depth due to a change in plan- 
etary radius between the line, Rpl, and the continuum, 
Rpc, can be related to the scale height, H, and the opac- 
ities in the absorption line, ki, and the continuum, n c ; 
this value can be approximated by the ratio of the opaci- 
ties multiplied by the area of an annulus one scale height 
thick relative to that of the stellar disk: 



{Rpl/R*) 2 — (Rpc/R*) 2 — 



2ttR p H 
ttR 2 



]n{Ki/ Kc ). (2) 



(jBrown et al. 1 12001[ ). To cause the observed transit 
depth difference in a hydrogen gas atmosphere would re- 
quire line opacity marginally greater than that of the con- 
tinuum (k;/k c ~8); the water-world composition would 
require an opacity that is unrealistically larger than that 
of the continuum (k;/k c ~2x 10 s ). Therefore from our 
CFHT/WIRCam observations, one would expect that 
the atmosphere of GJ 1214b must have a low mean 
molecular weight, a large scale height, and thus an at- 
mosphere dominated by hydrog en and/or helium. 

We compare our data to the [Miller- Ricci fc Fortnev I 
(2010) atmospheric models of GJ 1214b in Figure |3l 
We first employ the most conservative scenario and as- 
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Figure 6. Transit observations of the super-Earth GJ 1214b. The response functions are displayed inverted at the top of the plot 
at arbitrary scale. These are from left-to-right the MEarth (~0.78 fim), and Kitt-peak J-band (~1.25 fim) displayed in the red-curves, 
and the WIRCam J-band (~1.25 /ira), CH-40n filter (~1.69 fJ.m), and Ks-band (~2.15 fim) displayed in the black dotted curves. The 
Kitt-peak J-band bandpass overlaps the CFHT J-band bandpass. The VLT/FORS2 response functions are not included, as they essentially 
overlap with the points. The wei ghted mean of our J-band, Ks-band and ClrUOn filter observations are displayed by the black diamonds. 
The MEarth ( Charbonncau ct al. 2009), and the Kitt-peak Sada et al. (2010) transit depth a re displayed with the red diamo nds, while 
the VLT/FORS2 points IBean et al. I H2010T) are displayed with blue diamonds. We display two IMiller-Ricci fc Fortnev I 112010) GJ 1214b 
atmospheric models (a water- world model [grey dotted curve ] . and a solar- metallicity model [orange solid curve]); the other model we 
present in the cyan dashed curve is a Miller- Ricci & Fortncy (2010) solar-metallicity no methane model at wavelengths longer than one 
micron, while at shorter wavelengths we arbitrarily cut off the predicted absoprtion to simulate the impact of putative hazes. We integrate 
these atmospheric models over the WIRCam response functions and display these values in the appropriately coloured solid triangles. We 
caution that this comparison does not correct the depths for possible variability due to rotational modulation and thus any comparison 
should be considered illustrative, rather than definitive. 



sume that due to rotational modulation, or other sys- 
tematic errors, that we cann ot directly compare our mea - 
su red transit depths to the iCharbonneau et al~l (|2009l ) 
or IBean et al. I (|2010l ) depths. We thus scale-up or down 
the pr edicted absorption of the Millcr-R icci fe Fortnev I 
(2010) models by a multiplicative factor, analogous to 
an increase or decrease of the squared ratio of the plan- 
etary to stellar radius, to produce the best-fit (mini- 
mu m y 2 ) compared to our observ ations. We display 
the IMiller-Ricci fc Fortnev I (|2010t ) water-world model 
(an H2O dominated world) in the dotted grey line. 
We compare to two hydrogen-helium dominated atmo- 
spheres; the first has solar-metallicity (orange dot-dash 
line) while the second has solar metallicity but does 
not feature methane (cyan dot-dashed line). We inte- 
grate these models over the WIRCam response func- 
tions and calculate the associated x 2 °f our data com- 
pared to the model. Although, we cannot strongly dif- 
ferentiate between hydrogen-helium envelopes with so- 
lar metallicity (x 2 Q ; ar =9) from those without methane 
(xlo-methane= 7 )' the water-world composition is dis- 
favoured (Xi7 2 o = 26) by greater than 2a for our seven 
degrees of freedo m. Other high mean mol e cular weight 
models (e.g. the IMiller-Ricci fc Fortnev I (|20Tol ) C0 2 - 
dominated or H2 0/C02-mixture atmospheres) are dis- 
favoured with similar confidence. The support for the 
hydrogen/helium composition arises from the observed 
increased absorption in Ks-band as opposed to our J- 
band observations. The CB^On filter is nominally la 
discrepant from the observed models; it is unclear at this 
present time whether this discrepancy is physical or sim- 
ply due to low signal-to-noise. 

4.5. Comparison to observations at other wavelengths 

Arguably, completely excluding the constraints im- 
posed by transit depth observations at other wave- 
lengths, because of the effects of spots or systematic ef- 
fects, is unnecessarily conservative. Therefore, we also 



compare the weigh ted means of our transit depth mea- 
surements and the |Miller- Ricci fc Fortnev I (|2010l ) mod- 
els to the ICharbonneau et al. I (|2009r ) depth measure- 
m ent in the MEarth bandpass, the J-band measur ement 
of iSada et aT~l (|2010l) . and the IBean et"al~l (|2010D mea- 
surements from 0.78-1.0 /ira in Figure [6j We do not 
attempt to make a correction for the possibly variable 
spot activity between the various epochs at which the 
data were obtained; as discussed above in §4.21 rotational 
modulation can be expected to induce spurious transit 
depth variations as large as 0.03% of the ~1.35% transit 
depth of GJ 1214b near 1 \im, and will cause smaller vari- 
ations at longer wavelengths. A lso, a small discrep ancy 
will be induced by the fact that ISada et al. I (|2010P ) and 
ICharbonneau et al~~l (|2009l ) transit depth measurements 
were produced with the original ICharbonneau et al. I 

(2009) estimates of the inclination and other parame- 
ters fOTjJns_Sj^em i rather than the values derived from 
the IBean et al. I (|2010f ) white light photom etry that we 
use he re and that have been applied to the IBean et al. I 

(2010) spectrophotometry. For these reasons we caution 
that this comparison should be considered illustrative, 
rather than defini tive. 

Of the original IMiller-Ricci fc Fortnevl (|2010l ) atmo- 
spheric models the heavy mean molecular weight models, 
such as the water- world composition model (xlf 2 o = 40), 
are highly favoured over the low mean molecular weight 
compositions (Xs O ; ar =101 for the solar metallicity hy- 
drogen/helium dominated envelope) The water- world 
composition is favoured over the solar metallicity hydro- 
gen/helium dominated model with more than 5a confi- 
dence. Thi s is unsurpri s ing, a s it is largely the conclu- 
sion of the IBean et al. 1 (|2010f) paper, and results from 
the high precision of the VLT/FORS2 spectrophotom- 

14 As explained above in H4.4I we scale the radii of the planet in 
the models up or down to achieve the minimum \ 2 compared to 
the data. 
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etry and the lack of obs erved spectr a l feat ures in the 
very near- infrared. The iBean et al. 1 (|2010f ) paper ex- 
cludes the expected methane and water spectral absorp- 
tion features from a hydrogen/helium dominated atmo- 
sphere from 0.78 - 1.0 fim with high confidence. 

Another possibility to explain the la ck of observed 
spectral features are high altitude hazes (jFortnev 1120051) 
in the atmosphere of GJ 1214b that could mute the spec- 
tral features at shorter wave lengths. We discuss this pos- 
sibility further below ( §4.6j) . We thus also compare our 
observations to a no methane model where we cut off 
the absorption below 1 fim and set it equal to a nomi- 
nal value of 1.35%. At wavelengths greater than 1 fim, 
the values are identical to the no methane model. This 
abrupt transition is not intended to be physical, but sim- 
ply illustrative of the impact hazes could have on the ex- 
pected transmission spectrum at shorter wavelength;^. 
This model has the lowest x 2 (xlazy-h y drogen= M ) of an Y 
of the models as it is abl e to address t he lac k of observed 
spectral features in the IBean et al. I (|2010( ) spectropho- 
tometry, and our deeper Ks-band depth compared to our 
J-band deptrF^I. We note that the improvement in the \ 2 
of this haze model to the water- world model is not signif- 
icant, and both remain leading candidates to explain all 
the observations of GJ 1214 to date, as discussed below. 

4.6. Possible atmospheric compositions of GJ 1214b 

All the transmission spectroscopy observations of GJ 
1214b to date could be explained by a high mean molec- 
ular weight atmosphere if our deeper Ks-band transit is 
simply an outlier. If GJ 1214b's atmosphere does have a 
high mean molecular weight, a water vapour atmosphere 
is a leading possibility. This water- world scenario will re- 
main a viable candidate until observations are performed 
to either confirm that our Ks-band transit depths are in- 
deed deeper, or spectral features are detected at other 
wavelengths. 

A scenario that would qualitatively explain all the ob- 
servations to date, is a hydrogen/helium dominated at- 
mosphere with thick hazes that would mute the presence 
of spectral features arising from shorter wavelengths due 
to scattering. This haze layer would have to be at high 
altitudes, and low pressures (<200 mbar), to effectively 
mute the expected spectral features that would arise from 
being able to stare deep into the atmosphere of the planet 
in opacity windows in the very near-infrared. The ef- 
ficiency of scattering diminishes for wavelengths longer 
than the approximate particle size (|Hansen fc Travis"! 
119741 ) . The haze particles could not be much smaller than 
sub-micron size to acco unt for the the lack of observed 
spectral features in the IBean et al. I (|2010t ) spectropho- 

15 Depending on the size of the particles one would expect the 
observed transit radius to increase at very short wavelengths due 
to Rayleigh scattering or have a more complicated behaviour due 
to Mie scattering. 

16 While this p aper was in the late stages of revision, 
IDesert et al. I (|201 lr l presented Spitzer 3.6 and 4.5 /im channel ob- 
servations of the transit of GJ 1214b. They find similar tr ansit 
depths in these bands to t hose found by [Bean ct al. (20l3) and 
ICharbonneau et al. I (120091 ') in the very near-infrared. If methane 
is present in the atmosphere of GJ 1214b it should cause increased 
absorption in the 3.6 /im channel; as this is not observed, if the 
atmosphere of GJ 1214b is hydrogen/helium dominated, the com- 
bination of all the observations to date argues in favour of the hazy 
hydrogen/helium dominated model without methane. 



tometry in the very near-infrared. Due to the expected 
size of these putative haze particles, shorter wavelength 
optical observations would not be expected to show sim- 
ply the monotonic increase in planetary radius expected 
from a Rayleigh scattering signal, but instead the more 
complicated transmi ssion spectrum signal of Mie scatter - 
ing (see for example iLecavelier des Etangs et al~| [2008'). 

Such a haze or cloud layer is certainly not inconceiv- 
able. A cloud deck or haze has been reported to mute the 
optical transmission spect rum from 0.29 - 1.05 fim of the 
hot Jupiter HD 189733b (|Pont et al. |[200l [SnTFeTaTI 
\2QUh ; in the infrared HD 189733b may have absorption 
features with an opacity even greater than those that re- 
sult f rom the haze at those wavelengths (|Desert et al. I 
2010). The hazes of Jupiter and Titan may be other 
suitable analogies. Titan has a haze layer that is opti- 
cally thick in the optical, but has t ransparent windows as 
one moves into the near- infrared (|Tomasko et al. 1 ^008; 
Griffi tFet al. Ill993h - The opacity of hazes on Jupiter are 
high at short optical wavelengths, but are much smaller 
as one moves into longer optical wa velengths and into 
the near- infrared (jRages et al. II 1999h . 

A potential culprit for the particle causing this haze is 
a hydrocarbon derived from the photochemical destruc - 
tion of methane (|Moses et al. 120051: IZahnle et al. I2009D . 
Methane is found in the atmospheres of all the solar sys- 
tem's giant planets, as well as Titan. The end prod- 
uct of the breakdown of methane are higher order hy - 
drocarbons that condense as solids ([Rages et al. I fl9991 . 
Since GJ1214b should be relatively cool (T e ,-500K) its 
atmospheric carbon inventory could feature abundant 
methane, like Jupiter. Particulates with a relatively 
small mixing ratio can have important effec ts at the 
slant viewing geometry appropriate for transits (jFortnev I 
120051 ) , and thus hydrocarbons in a high altitude haze are 
one possible exp l anatio n the lack of observed features in 
the IBean et al. I (|2010D spectrophotometry. 

We lastly note, that the actual spectral features of GJ 
1214b, whether its atmosphere is hydrogen/helium dom- 
inated or not, coul d be more complicated and thu s very 
different than the Hiller-Ricci fe Fortnevl (|2010f ) mod- 
els predict. One such reason could be the impact of 
non-equilibrium chemistry, which will be explored in a 
forthcoming publication (Miller-Ricci Kempton et al. in 
prep.). 

Clearly further observations are required to differenti- 
ate between these scenarios and determine the true at- 
mospheric makeup of GJ 1214b. 

4.6.1. An opacity source at ^2.15 fim 

The increased transit depths we note in our Ks-band 
observations, argue for an opacity source near —2.15 
microns that is causing increased absorption along the 
limb of the planet. One such, possible opacity source 
is methane, which is pre dicted to cause absorption from 
-2.2 to -2.4 fim in the iMiller-Ricci fc Fortnevl ()2010f ) 
hydrogen-helium dominated atmospheric models. Al- 
though our Ks-band transit depth is qualitatively con- 
sistent with this methane spectral absorpt i on fea ture, 
we note that the IMiller-Ricci fc Fortnevl (|2010l ) hy- 
drogen/helium model with solar metallicity but with- 
out methane, provides a near-identical goodness-of-fit as 
compared to the solar metallicity model with methane. 
This is because the methane absorption feature is present 
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at the red-edge of our Ks-band only; on the blue edge of 
the Ks-band the hydrogen/helium model with methane 
actually features less absorption than the model with- 
out methane, so overall the predicted transit depth is 
similar whether methane is or is not present according 
to the iMiller-Ricci k Fortnev"! (|201Q[ ) prediction. Also, 
greater concentrations of methane in the atmosphere are 
not expected to cause increased absorption at these wave- 
lengths. In the no methane model, the increased Ks- 
band absorption compared to J-band results from wa- 
ter opacity. Thus both water and methane remain vi- 
able candidates for this increased absorption. Both these 
molecules should also lead to spectral featu res from 0.78 
- 1.0 J im that have been ruled out by the IBean et al. I 
(2010) spectrophotometry at high confidence. Thus for 
these chemicals to remain viable opacity sources, one re- 
quires the presence of hazes in the atmosphere of GJ 
1214b, o r that the spectral feature s are m ore complicated 
than the IMiller-Ricci k Fortnevl (|2010D models predict. 
We encourage further modelling to elucidate whether this 
^2.15 \xm feature is due to methane, water or another 
opacity source. 

4.7. Consequences of a hydrogen/helium dominated 
atmosphere 

A hydrogen/helium dominated atmosphere on GJ 
1214b would be expected to undergo significant hydrody- 
namic escape. Therefore, if GJ 1214b is hydrogen/helium 
dominated then it may have lost or is losing a significant 
fractio n of its gaseous enve l ope. ICharbonneau et al~l 
(|2009fl and lRogers fc Seager I (|2010T > predicted that if GJ 
1214b's atmosphere is dominated by hydrogen gas then 
it will lose on the order of ~10 9 gs" 1 , or ^0.02 M ffi on a 
4 Gyr timescalj 17 !. As its host star may have been more 
active earlier in its life, and thus brighter in the ultravio- 
let, its cumulative mass loss may be higher. Thus, either 
if GJ 1214b's hydrogen/helium envelope is primordial or 
due to outgassing we may expect it to have lost and will 
be losin g a non-neglig i ble fra ction of its atmosphere. 

Also, iCarter et al. I (|2011[) have pointed out that the 



fit the observ ed mass and radius - and thus density 
constraints of ICharbonneau et al. I (120091). We com pare 



radius of the star, GJ 1214, that one obtains from stel- 
lar evolutionary models is very different than the radius 
one obtains from an analysis of the light curve. As a 
result the density of GJ 1214b varies accordingly, from 
one where a significant gaseous atmosphere is likely, to a 
much higher density where one would only expect a thin 
gaseous atmosphere on top of a solid terrestrial planet. 
Our increased Ks-band transit depth argues in favour of 
the lighter density and the stellar radius suggested by fits 
to light curve parameters. 

4.7.1. Constraints on GJ 1214b's bulk composition 

If the atmosphere of GJ 1214b is hydrogen/helium 
dominated then this allows us to place a constraint on 
the planet's bulk composition, namely its core, man- 
tle and possibly its ice layer ((Rogers k Seager 1 120101 : 
iNettelmann et al~l 120101) . This is because the lighter 
atmospheric composition of a hydrogen/helium atmo- 
sphere, compared to for instance a water-world com- 
position, requires a heavier interior composition of sil- 
icates, ferrous material or ices to compensate in order to 



17 ICh arbonneau ~et al. I l|2003 ') report a 3-10 Gyr age for this 
system 



our results to two numerical m o dels ((Rogers k Seager I 
120101 INettelmann et al. I l2010h that attempt to deter- 
mine the range of realistic bulk compositions that agree 
with the observed mass and radius constraints. Regard- 
less of whether the planet's hydrogen/helium envelope 
is primordial or due to outgassing, it is expected that 
this atmospheric layer will be a small percentage ( ~5%) 
of this planet's t otal mas s dRogers k Seager I 120101 : 
INettelmann et al. I 12010ft. INettelmann et al. I (|2010l ) 
and lRogers k Seager I (|20iol ) predict, under their hydro- 
gen/helium atmosphere scenarios, that a wide range of 
core/mantle masses is still viable (from a few percent to 
~99% of the mass of GJ 1214b). The higher core/mantle 
masses result from a planet with very little water, while 
for the lower m asses it would enta il a ma ssive interior 
water /ice layer. INettelmann et al~l ()2010h suggest that 
if this planet's atmosphere is dominated by hydrogen and 
helium then one can place an upper limit on the water to 
rock ratio of approximately six-to-one; the true value of 
this quantity is expected to be much lower, and thus the 
core mass is expected to make-up at least ~14% of the 
plane t's mass and likely much more (iN ettelmann et aT~l 

I2010I) . 

5. CONCLUSIONS 

We report observations of four transits of GJ 1214b 
using WIRCam on CFHT. We observed nearly simul- 
taneously in J-band and in Ks-band for three of the 
transits, and in J-band and the CILiOn filter in an- 
other. Our best-fit J-band transit depth is consis- 
tent with the values obtained in the optical and very 
near-infrared: (i? P .//.R*) 2 = 1 -338±0.013%. Our Ks-band 
transit is deeper: (R PKs /R t ,) 2 =lA38±0M9%. Our J 
and Ks-band transit depths are inconsistent at the 5cr 
level. The impact of spots and rotational modulation 
on the transit depths we observe should be small; nev- 
ertheless spots will induce small changes in the transit 
depths we measure from epoch to epoch, and as a re- 
sult a better metric to quantify our observations may be 
the factor that the Ks-band transits are deeper than the 
J-band transits observed simultaneously. Our Ks-band 
transits display a deeper depth than our J-band transits 
by a factor of \r PKs /Rpj) 2 =1.072±0M8. We thus de- 
tect increased transit depths in our broadband Ks-band 
(^2.15 /um) as compared to J-band (~1.25 /an) with 4cr 
confidence. The difference in transit depth between the 
two bands that we measure is best explained due to a 
spectral absorption feature from the atmosphere of GJ 
1214b; the only way to get a spectral absorption fea- 
ture this prominent is if the atmosphere of GJ 1214b 
has a large scale height, low mean molecular weight and 
is thus hydrogen/helium dominated. Water or methane 
are possible opacity sources to explain this absorption. If 
GJ 1214b's atmosphere is hydrogen/helium dominated a 
range of core/mantle masses and ice layers is still viable, 
but the planet must have a large rocky core/mantle in- 
terior to its gaseous envelope. In this case, our increased 
Ks-band transit depth would be the first detection of a 
spectral feature in a super-Earth atmosphere, and GJ 
1214b would best be described as a mini-Neptune. 

However, when combining our observations with other 
observations of GJ 1214b, most specifically the lack 
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of spectral features observed in the iBean et al. I (|2010l ) 
VLT/FORS2 spectrophotometry from 0.78 - 1.0 /im, the 
atmospheric composition of GJ 1214b is less clear. There 
are several leading possibilities. One possibility remains 
that the atmosphere of GJ 1214b is water-vapour dom- 
inated and our increased Ks-band transit depth is sim- 
ply an outlier; our increased Ks-band transit depth will 
have to be reconfirmed or spectral features at other wave- 
lengths will have to be detected before this scenario can 
be confidently ruled out. The possibility that is arguably 
the most consistent with all the observed data to date, 
is that GJ 1214b has a hydrogen/helium dominated at- 
mosphere with a haze layer at high altitude consisting 
of particles that can be no smaller than approximately 
sub-micron in size; such a scenario would explain the 
lack of observed sp ectral absorp t ion fe atures in the very 
near-infrared in the IBean et al. I (|2010l ) spectrophotome- 
try and our own increased Ks-band transit depth if there 
is an opacity source at ~2.15 fmi. Lastly, we note that 
the true spectrum of GJ 1214b could be more compli- 
cated than our models predict for a variety of reasons - 
one such possibility is the importance of non-equilibrium 
chemistry in GJ 1214b's atmosphere which would al- 
ter GJ 1214b's predicted transmission spectrum (Millcr- 
Ricci Kempton et al. in prep.). 

Clearly, further observations are required to precisely 
determine the nature of GJ 1214b's atmosphere. We en- 
courage further observations to confirm our increased Ks- 
band transit depth. We plan to use CFHT/WIRCam to 
reobserve the transit of GJ 1214b in Ks-band on several 
occasions in the Spring/Summer 2011 observing season. 
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